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Abstract 


In the framework of heavy quark effective theory we use QCD sum rules to 
calculate the masses of the és (0*,1*) and (1,27) excited states. The results are 
consistent with that the states D,7(2317) and D,7(2460) observed by BABAR and 
CLEO are the 0T and 1* states in the jı = Iy doublet. 


PACS number: 12.39.Hg, 13.25.Hw, 13.25.Ft, 12.38.Lg 


1 Introduction 


Recently BaBar Collaboration announced a positive-parity narrow state with a rather 
low mass 2317 + 3 MeV in the Dł}(1969)r channel [1], which was confirmed by CLEO 
[2] and BELLE [3] later. Because of its low mass and decay angular distribution, its 
J? is believed to be 0*. In the same experiment CLEO [2] observed a state at 2460 
MeV with the possible spin-parity J? = 1+ in the D*m channel. BaBar [1] also found a 
signal near there. Since these two states lie below DK and D*K threshold respectively, 
the potentially dominant s-wave decay modes D,7;(2317) —> D,K etc are kinematically 
forbidden. Thus the radiative decays and isospin-violating strong decays become favorable 
decay modes. The later decay goes in two steps with the help of virtual Dsn intermediate 
states, D,;(2317) — Dsn — Dsn? where the second step arises from the tiny 7-7° isospin- 
violating mixing due to m, 4 mq. 

The experimental discovery of these two states has triggered heated debate on their 
nature recently in literature. The key point is to understand their low masses. The 
D,,7(2317) mass is significantly lower than the values of 0* mass in the range of 2.4 — 2.6 
GeV calculated in quark models of [4]. The model using the heavy quark mass expansion 
of the relativistic Bethe-Salpeter equation in [5] predicted a lower value 2.369 GeV of 0+ 
mass which is still 50 MeV higher than the experimental data. Bardeen, Eichten and 
Hill interpreted them as the és (0*,1*) spin doublet as the parity conjugate states of the 


(07,17) doublet in the framework of chiral symmetry|6]' (See also Ref. [7]). A quark- 
antiquark picture was also advocated by Colangelo and Fazio [9], Cahn and Jackson [10], 
Godfrey [11]. Based on such a ”conventional” picture, the various decay modes were 
discussed in Refs. [6, 9, 11]. 

Apart from the quark-antiquark interpretation, D, ;(2317) was suggested to be a four- 
quark state by Cheng and Hou [12] and Barnes, Close and Lipkin [13]. Szczepaniak even 
argued that it could be a strong Dz atom [14]! But we think it would be very exceptional 
for a molecule or atom to have a binding energy as large as 40 MeV. 

Van Beveren and Rupp [15] argued from the experience with ag/ fo(980) that the low 
mass of D,7(2317) could arise from the mixing between DK continuum and the lowest 
scalar nonet. In this way the 0* és state is artificially pushed much lower than that 
expected from quark models. 

The recent lattice calculation suggests a value around 2.57 GeV for 0* state mass 
[16], much larger than the experimentally observed D, (2317) and compatible with quark 
model predictions. The conclusion in Ref. [16] is that D,;(2317) might receive a large 
component of DK and the physics might resemble ag/ fo(980). Such a large DK compo- 
nent makes lattice simulation very difficult. 

In this paper we shall use QCD sum rules [20] in the framework of the heavy quark 
effective theory (HQET) [21] to extract the masses since HQET provides a systematic 
method to compute the properties of heavy hadrons containing a single heavy quark via 
the 1/mg expansion, where mg is the heavy quark mass. The masses of ground state heavy 
mesons have been studied with QCD sum rules in HQET in [22, 23, 24]. In [17, 18, 19] 
masses of the lowest excited non-strange heavy meson doublets (0*,1*) and (1*, 2+) were 
studied with the sum rules in HQET up to the order of O(1/mg). These masses were 
also analyzed in the earlier works [25] with sum rules in full QCD . In this short note we 
extend the same formalism in [17, 18] to include the light quark mass in calculating the 
cs mesons. We shall also use more stringent criteria for the stability windows of the sum 
rules in the numerical analysis. 


2 my—mg at the leading order of HQET 


The proper interpolating current J ba for a state with the quantum numbers j, P, je in 


HQET was given in [17]. These currents were proved to satisfy the following conditions 


(Oleg. (OI P's je) = Sersjydpp pgn , (1) 
j a lg a j a1 ajbi 
1 (OJT (E (x) Jj py (0)) |0) = jj OP P10 5,31(—1)? Sg; Ar EET 
x / dt5(a — vt) Tp,,(2) (2) 


in the mg — oo limit, where 7%" is the polarization tensor for the spin j state, v 
is the velocity of the heavy quark, gẹ? = g°® — v%v® is the transverse metric tensor, S 


'The existence of parity doublets has also been shown by combining chiral symmetry and heavy quark 
symmetry in Bethe-Salpeter approach in [8]. 


denotes symmetrizing the indices and subtracting the trace terms separately in the sets 
(a, ---a,;) and (G;---G;), fej, and IIp;, are a constant and a function of x respectively 
which depend only on P and jg. 

We consider the correlator 


: ika aya; Brb 
Mw) =i f atze"? (0T (155 a (0)) 10) (3) 
where w = 2u-k. It can be written as 
a fie x 
Hw) = ———+—_ + higher states , (4) 
2A; Pa, = 0 


Where Aj Pje = liMmo—>o(Mm; pz — Mo). On the other hand, it will be calculated in 
terms of quarks and gluons. Invoking Borel transformation to Eq. (3) we get 
2A j, P,j 1 We w 
fjue a p(w)e T + condensates . (5) 
Á T J2mMmq 

where m, is the light quark mass, p(w) is the perturbative spectral density and we is 
the threshold parameter used to subtract the higher state contribution with the help of 
quark-hadron duality assumption. 

We shall study the low-lying (0+, 1%), (1+, 2+) (és) states. The values of various QCD 
condensates are 


(gs) = —(0.8 + 0.1) * (0.24 GeV)? , 
(a,GG) = 0.038 GeV* , 
mê = 0.8 GeV?. (6) 


We use m,(1GeV) = 0.15 GeV for the strange quark mass in the MS scheme. We use 
Agcp = 375 MeV for three active flavors and Agcp = 220 MeV for four active flavors. 
The sum rules with massless light quarks have been obtained in [18, 17]. 

In the numerical analysis of the QCD sum rules we require that the high-order power 
corrections be less than 30% of the perturbation term. This condition yields the minimum 
value Tmin of the allowed Borel parameter. We also require that the pole term, which is 
equal to the sum of the cut-off perturbative term and the condensation terms, is larger 
than 60% of the perturbative term, which leads to the maximum value Tnax of the allowed 
T. Thus we have the working interval Tmin < T < Tmax for a fixed we. If Tmin > Tmar, 
we are unable to extract useful information from such a sum rule. In the ideal case, 
the difference between the meson mass m jy and the heavy quark mass mg (or other 
observables) is almost independent on T for certain values of we. Namely, the dependence 
on both the Borel parameter and continuum threshold is minimum. In realistic cases, the 
variation of a sum rule with both T and w, will contribute to the errors of the extracted 
value, together with the truncation of the operator product expansion and the uncertainty 
of vacuum condensate values. 

In subsections 2.1 and 2.2 we give sum rules for (0*,1*) and (1+, 2*) doublets using 
interpolating currents with derivatives, respectively. For comparison we discuss sum rules 
for (0*,1*) using interpolating currents without derivatives in subsection 2.3. 
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Figure 1: The variation of A(}") (in unit of GeV) of the (0T,1*) doublet with T and 
We for the derivative currents. The vertical and horizontal axes correspond to A and T. 
From top to bottom, the curves correspond to we being 3.1, 2.9, 2.7 GeV respectively. 


2.1 (0°,17) Doublet With Derivative Currents 
We employ the following interpolating currents [17] 


, T 2 À 
Tics = AT h,(—i) Pis , 


_ lS eel 
jhe = Z hoye (i) Pes - (7) 


where D¥ = D! — (v - D) v”, with D” = 0" — ig A” being the gauge-covariant derivative. 


We obtain the sum rule relevant to A 


= 3 We 
fae E f lwt + 2m? — 6m2w? — 12mwe™ T dw 
1 B 3 7 Ms 
= me (5s) + gs (8s) — ic: (as G’) . (8) 


Taking the derivative of the logarithm of the above equation with respect to (£) one 
obtains the sum rule for A. Substituting the obtained value of A in (8) one obtains the sum 
rule for f’, 1/2 In Figure 1, the variation of A(1/2*) with T and w, is shown. According 
to the criteria stated above the working range is 0.38 < T < 0.58 GeV. We have 


A(1/2*) = (0.86+0.10) GeV, 
fia = (0.314 0.05) Gev?”, (9) 


where the central value corresponds to T = 0.52 GeV and we = 2.9 GeV . 
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Figure 2: The variation of Aa) (in unit of GeV) of the (17, 2+) doublet with T and we. 
From top to bottom, we = 3.2, 3.0, 2.8 GeV respectively. 


2.2 (1,2) Doublet 
For the (1*,2*) doublet, by using the following interpolating currents [17] 


a 35 A Q 1 a 
ia = [Erri (2; = 3T P) q, 


a1,a 1, (-1 a a a a 2 aja 
Bai > 2 hy | 9 ) Gz. Pi = 3 Ho” Pr) q, (10) 
we obtain the sum rule 
= 1 We 
F gae R = ao | [wt + 2m,w? — mw? — 12mĉwle T dw 
i T 2Ms 
l1 3 1 & 1 Ms ,O 
= mo (88) — O T Sa . 11 


¿From the sum rule, the variation of A(3/2*) with T and we is plotted in Figure 2, . 
The working range is 0.55 < T < 0.65 GeV. We have 


A(3/2*) = (0.83 + 0.10) GeV, 
frre = (0.194 0.03) GeV*?, (12) 


where the central value corresponds to T = 0.62 GeV and w, = 3.0 GeV . 
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Figure 3: The variation of the ratio of the pole contribution and the perturbative piece of 
the (0*,1*) doublet with T and we for the non-derivative currents. From top to bottom, 
We = 3.1, 2.9, 2.7 GeV respectively. 


2.3. (0+,1*) Doublet With Currents Without Derivative 


Possible different currents for (0*,1*) doublet are the non-derivative currents 


With the non-derivative currents the sum rule reads 


_on 3 we o 
Tane NT a: TE ~ W? — 2m,w — 2m?je™ T dw 
1 5 1 
+3 (Ss) + z (5s) — 372 me (3s) . (14) 


Requiring that the condensate contribution is less than 30% of the perturbative term, 
we get Tmin = 0.75 GeV. In Figure 3, we show the ratio of the pole term, i.e the sum of 
the cut-off perturbative term and condensation terms, for the sum rule (14). In the whole 
range of T > Tmin, this pole contribution is less than 40%. Hence, There is no stability 
window for this sum rule satisfying our criteria stated before. 

If we arbitrarily loosen the analysis criteria and require the pole contribution is greater 
than 30% only, we get the working range 0.75 < T < 1.2 GeV. In Figure 4, the variation 
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Figure 4: The variation of A (in unit of GeV) of the (0+,1*) doublet with T and we for 
the non-derivative currents. From top to bottom, we = 3.1, 2.9, 2.7 GeV respectively. 


of A(1/2*) with T and w, is shown. Numerically, 


A(1/2*) = (1.30 + 0.15) GeV, 
fraye = (0.39 + 0.05) Gev*?, (15) 


where the central value corresponds to T = 1.0 GeV and we = 2.9 GeV. Because of 
the weaker criteria used these results are less reliable and will not be used in the final 
numerical results. 

Here we would like to make some remarks. Usually, the currents with the least number 
of derivatives are used in QCD sum rule approaches. The sum rules then are less sensitive 
to the threshold energy we. However, it is pointed out in [17] that in the non-relativistic 
limit the coupling constant of these currents to the P wave states vanishes. If this coupling 
constant is suppressed due to this reason, the relative importance of the contribution of 
the DK and other states in continuum in the sum rules which are usually neglected 
would be enhanced. Besides, it is shown in [26, 27] by using soft pion theorem that the 
contribution of the Dr continuum is large in the sum rule with the non-derivative current 
for the 0+ state of the non-strange D system and significantly decreases the value of A. 
Similar method of calculation is not good in the present case, but it indicates that the 
contribution of the DK continuum with the non-derivative current may be large too. 


3 Sum rules at the 1/mg order 


To the order of 1/mg, the Lagrangian of HQET is 
na 1 1 P 
Leg = h iv- D hy +—K + —S+ O(1/m9), (16) 
2MQ 2MQ 
where h,(a) is the velocity-dependent effective field related to the original heavy-quark 
field Q(x) by 
. 1 
h(a) = emare LEE Qn). (17) 

K is the kinetic operator defined as 


S is the chromo-magnetic operator 


$= 5 CmaglMQ Ju) Pay Opp GP” Py , (19) 
3/Bo 
a 
See re 
ans) i 


Considering O(1/mg) corrections the pole term of the correlator on the hadron side 
becomes 


where Cmag = ( 


(f tof)" ie 25m f’ 2fof 
I öle “AYA a aN AR Fok ND ak 2 
(opa Arma Ase a Whew 
where ôm and ôf are of order O(1/mg). 
To extract dm in (20) we follow the approach of [23] to consider the three-point corre- 


lation functions 


J J 


dolia (ww!) = 2 f Peagi [0) , (21) 


where O = K or S. The scalar function corresponding to (21) can be represented as the 
double dispersion integral 


1 d 1 
olw, w) = = / nee es 


nm J (s—w)(s'—w!) ` 
The pole parts of doIH(w,w’),O = K, S are 


PK Pado) POl) 


Ul ! ee = peel aes 
KTW, "poe QA—w)2A-—w') 2A-—w 2A —u! ’ 2) 
fe 4 ke duf’ X 2/Gs(w') | Gs(w) 
eee OK ahaa) |haw oh) 


where 

K;,pi = (J, P, july (D1)" hold, P, ji) 5 

2dM¥j Pj = (j, P, jiho Jo, G” halj, P, ji) , 

dm = dijs Gj—rjeg = 2j +2, dj41/2j = 2j 
Letting w = w’ in eqs. (22) and (23) and comparing with (20), one obtains [23] 

1 
ôm = ———(K + dmCmag®) . (24) 
4mo 


The single pole terms in (22) and (23) come from the region in which s(s') = 2A and 
s'(s) is at the pole for a radial excited state or in the continuum. They are suppressed 
by making the double Borel transformation for both w and w’. The Borel parameters 
corresponding to w and w’ are taken to be equal. One obtains thus the sum rules for K 
and X as 


PK pram — f l l dwdw e i IPT on (w, w") ’ (25) 
2Ms J 2Ms 
PENT f] f dwd tT psw, u’). (26) 
2Ms J 2Ms 
In this section we shall neglect the m, corrections to K and X. We obtain for the 
7 = 5 doublet with derivative currents 
e 3 f“ = 3 
fe aK e a = Fol we a tae 
1 
De PAT — a f o e-+ dsds! + —(a,GG 27 
pk ae sds! + = —(0,GG) T (27) 
with 


1 a,(2T) 


OS 


We can eliminate the dependence of the K and X on f4 1/2 and A through dividing the 
above sum rules (27) by the sum rule (8). We use the same working windows as those of 
two-point sum rules for A in evaluating O(1/mg) corrections. Numerically we have 

Ky/2 = (—1.60 + 0.30) GeV? 
Y1/2 = (0.28 + 0.05) GeV?. (28) 


1+ 
The variations of K and © with T and we for the j? = 5 doublet are shown in Figure 5 
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3 
For jf = 5 doublet, we have 


a 1 fue 5 _, 7 
Fopke 2A/T _ 9778 o wê e™®/T dw + TE T3 
1 
2 —2Ā/T _ — 
Pane =f f o “SF dsds! + —— ——(asGG) T (29) 
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Figure 5: The variation of Kı;2 (in unit of GeV’) of the (0+,1*) doublet with T and we 
for the derivative currents. From top to bottom, we = 3.1, 2.9, 2.7 GeV respectively. 
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Figure 6: The variation of X1/2 (in unit of GeV?) of the (0*, 1+) doublet with T and we 
for the derivative currents. From top to bottom, we = 3.1, 2.9, 2.7 GeV respectively. 
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Figure 7: The variation of K3/2 (in unit of GeV”) of the (1+, 2+) doublet with T and we. 
From top to bottom, we = 3.2, 3.0, 2.8 GeV respectively. 


with 


Cu s= (s = 3s') + (s = s')}. 


at 
The variations of K and E with T and we for the jë = 5 doublet are plotted in Figure 


7 and 8 respectively. Numerically we have 


K3)2 = (—1.64 + 0.40) GeV? 
Y3/2 = (0.058 + 0.01) GeV? . (30) 


The spin-symmetry violating term S not only causes splitting of masses within the 
same doublet, but also causes mixing of states with the same j, P but different jı. In Ref. 
[18] corrections from the mixing is found to be negligible. We omit this effect here. 


4 Numerical Results and Discussions 


We present our results for the cs system assuming the HQET is good enough for excited 
D mesons. For the doublet (0+,1*) with the derivative current, the weighted average 
mass is 


1 1 
= (mp. + 3mp» ) = me + (0.86 + 0.10) + — [(0.40 + 0.08) GeV’). 
4 s0 a Me 
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Figure 8: The variation of 3/2 (in unit of GeV”) of the (1*,2*) doublet with T and we. 
From top to bottom, we = 3.2,3.0, 2.8 GeV respectively. 


The mass splitting is 


1 
mp: — mp = — [(0.28 + 0.05) GeV’) . 


Cc 


For the (1+, 2+) doublet we have 


1 1 
5 (3mp,, + 5mp=,) = Me + (0.83 + 0.10) + aa [(0.41 + 0.10) GeV?) . 


(31) 
The 1*,2* mass splitting is 


1 
Mp*, — Mp, = — [(0.116 + 0.06) GeV] . 
s Me 
The results for bs system are obtained by replacing me by m, and multiplying © by 
0.8 (since Cmag © 0.8 for B system by using the values of Agcp given in Sect. 2) in above 
equations. 
Choosing m, to fit the experimental value 


1 
= (8mp,, + 5Mp«,) = 2.56 GeV, 
8 sl s2 


where we again neglect the mixing between two 17 states, we obtain me = 1.44 GeV. 
Using this m, value we obtain the following numerical results. The 1*,2* mass splitting 
is 


mp:, — Mp, = (0.080 + 0.042) GeV 


12 


which is consistent with the experimental value 37 MeV within the large theoretical 
uncertainty. Experimentally, the mass splitting in the (17,2+) doublet in D, system 
is almost equal to that in D system. This justifies neglecting the mMm, correction to the X 
term in our calculation. For (0*,1*) doublet, the weighted average mass is 


1 


The mass splitting is 
Mp*, — mps, = (0.19 + 0.04) GeV . (32) 


Therefore, the 0+ mass is predicted to be mo+ = (2.42 + 0.13) GeV. This is consistent 
with the experimental value 2.317 GeV, though the central value is 100 MeV larger than 
data. If the 1* D, state of mass 2.460 GeV found by CLEO [2, 1] is assigned as the other 
member of the (0*,1*) doublet, then the observed mass splitting 0.143 GeV is consistent 
with (32). Notice that the 41/2 value in Fig. 6 changes very slowly between T = 0.35—0.8 
GeV. Therefore the result (32) is insensitive to the stability window used. 

We would like to note that the stability of the sum rules (in particular, those for the 
K and X) obtained by us is not as good as those for the ground states, as can be seen 
from the figures. And for the charm flavor hadrons, the 1/m, corrections are significant. 
So the predictions on masses have a relatively large uncertainties, which are estimated by 
given errors. 

If we replace the strange quark condensate by up/down quark condensate and m, by 
the zero up/down quark mass, we can extract the non-strange excited D meson masses. 
For (1+,2*) doublet, the sum rule window is 0.57 < T < 0.67 GeV for we = 2.8 — 3.2 
GeV. The results are 


K37. = —(1.6 + 0.4) GeV’, 
Y3/2 = (0.057 + 0.01)GeV” (33) 


with the central value at T = 0.62 GeV. For (0*,1*) doublet, the sum rule window is 
0.46 < T < 0.6 GeV for we = 2.7 — 3.1 GeV. The results are 


Ky. = —(1.57 + 0.4) GeV’, 
Y1/2 = (0.286 + 0.05)GeV? (34) 


with the central value at T = 0.53 GeV. The effect of strange quark is to make the mass 
of D, mesons a little bit larger, as expected. Experimentally, Belle collaboration found 
Mp, = (2290 + 22 + 20) MeV and mp; = (2400 + 30 + 20) MeV recently [29] while 
CLEO collaboration found mp: = 2400745 MeV [30]. The j3/2 non-strange doublet mass 
is known precisely [28], mp, = 2460 MeV and mp, = 2420 MeV. Our results (33)-(34) 
are consistent with the experimental data within theoretical uncertainty. 

In summary, we have calculated the masses of the excited (0*,1*) and (1*,2+) dou- 
blets for the és system to the 1/mg order in the HQET sum rules. The numerical results 
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imply that the D,;(2317) and D,7(2460) observed by BABAR and CLEO can be consis- 
tently identified as the (0*,1*) doublet with jı = within the theoretical uncertainty. 
Especially, the mass splitting in the (0*,1*) doublet is reproduced quite well. The re- 
pulsion between these states and the Dk and DK* continuum may help to lower their 
masses. In the framework of the sum rule this effect should come from the contribution 


from the Dk and DK* continuum to the dispersion integral. 
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